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ABSTRACT 

This is the second in a series of papers presenting an analytical model for the evolu- 
tion of FRII radio sources. In this paper we evaluate the expected radio emission from 
a radio source incorporating energy loss processes for the relativistic electrons. By 
combining these results with our earlier dynamical model we calculate evolutionary 
tracks through the Power-Linear size diagram. These tracks are in good agreement 
with the observed distribution of sources in this diagram. The effects of different forms 
for the evolution of the magnetic field in the cocoon, the redshift of the source, the 
environment of the source and the defining parameters of the jet are investigated. 
The evolutionary tracks are found to be insensitive to the assumed form of the mag- 
netic field evolution. Some evidence against protons as a major constituent of the jet 
material is also found. 

Key words: Radiation mechanisms: cyclotron and synchrotron - Galaxies: active - 
Galaxies: jets 



1 INTRODUCTION 

This is the second in a series of papers presenting an analyti- 
cal model for the evolution of FRII radio sources; in the first 
paper (Kaiser & Alexander 1997, hereafter KA97) we inves- 
tigated the dynamical evolution of the source and here we 
extend this analysis to consider the luminosity evolution of 
a radio source and the tracks it follows through the Power- 
Linear size (P-D) diagram. The P-D diagram was intro- 
duced by Shklovskii (1963) as a powerful tool for investigat- 
ing the temporal evolution of FRI and FRII radio sources. 
By plotting the radio luminosity at a specific frequency, P v , 
as a function of the linear size of a source, D, a diagram 
analogous to the Hertzsprung-Russel diagram is obtained. 
The evolution of a radio source through the P-D diagram 
is, however, not well understood. Some constraints can be 
placed on these evolutionary tracks by the relative densi- 
ties of objects in regions of the P-D diagram; specifically, 
sources with large linear sizes (D > 1 Mpc) and high radio 
luminosities (P„ > 10 26 W Hz" 1 sr" 1 at v = 178 MHz) seem 
to be rather rare suggesting that the luminosity of sources 
should decrease quickly with linear size for sizes approaching 
1 Mpc. Care must be taken in interpreting the P-D diagram 
in this way since selection effects have a strong influence on 
the observed high-redshift population (e.g. Masson 1980, 
Macklin 1982) and giant sources are likely to be difficult 
to detect since the extended emission of the cocoon will be 
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weak and the hot spots may therefore appear unconnected 
(e.g. Muxlow & Garrington 1991). Searches specifically for 
giants have been undertaken by Cotter et al. (1996) and 
Subrahmanyan et al. (1996). Both groups identify a number 
of radio sources in low frequency surveys as classical doubles 
of extreme size, although almost all of these sources have lu- 
minosities below 2 x 10 26 W Hz" 1 sr -1 suggesting that there 
is a real deficit of powerful giant sources. 

Baldwin (1982) used the model for FRII sources of 
Scheuer (1974) and minimum energy arguments (e.g. Bur- 
bidge 1956) to calculate evolutionary tracks through the P- 
D diagram. In its original form his model predicts P oc D 7 / 8 
for a power law energy spectrum of the relativistic elec- 
trons, N(E)dE oc E~ p dE, with p = 2 as exponent. This 
implies that all sources would crowd together in exactly 
that area of the diagram which is almost empty. Modifi- 
cations of the model incorporating a density gradient in the 
external medium surrounding the radio source (p x oc d _/3 , 
where d is the distance from the center of the distribution) 
predicted tracks with decreasing radio luminosity in reason- 
able agreement with observations for /3 = 1.9 and D < 300 
kpc. Lo explain the rapid decline in luminosity for D ~ 1 
Mpc Baldwin suggested that the density of the atmosphere 
may decrease more steeply with /3 — 2.9 for D > 300 kpc. 
However, sources in environments with a density gradient 
steeper than 1/d 2 are unlikely to develop the typical FRII 
morphology (Falle 1991, Kaiser & Alexander 1997). 

An alternative hypothesis to account for the decrease 
in luminosity at large D is that the central engine ceases 
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to supply energy to the jet after a certain period which is 
similar for all sources (for a discussion see Daly 1994). The 
cocoon will then dim quickly via adiabatic expansion. Ob- 
servations indicate that only of order 4% of all radio sources 
show no signs of ongoing AGN activity Giovannini (1988) in 
support of this hypothesis. 

The synchrotron radiation of FRIIs demonstrates the 
presence of highly relativistic electrons in the cocoon (e.g. 
Rees 1971). To achieve charge neutrality in the cocoon we 
need some species of positively charged particles. If we as- 
sume that there is no entrainment across the contact discon- 
tinuity which is separating the cocoon from the IGM the jet 
must supply these additional particles. The nature of these 
charge-balancing particles is subject to debate (e.g. Leahy 
1991). If they are protons, a significant fraction of the kinetic 
energy transported by the jet will be stored by these in the 
cocoon but will not contribute to the radiated power. Al- 
ternatively, if the particles providing charge neutrality were 
positrons, they would behave almost exactly like electrons 
in a completely tangled magnetic field and will contribute to 
the observed synchrotron emission. The kinetic power of the 
jet required to produce a cocoon of given radio luminosity 
would be significantly reduced in the case of a pair plasma. 

In this paper we extend the model described in KA97 
for FRII sources by the addition of energy loss processes of 
the relativistic electrons in the cocoon to compute evolution- 
ary tracks through the P-D diagram. We show that changes 
in the external atmosphere or a mechanism for switching 
off the central engine are not required to reproduce the ob- 
served distribution of sources in the P-D diagram. We also 
find some evidence against protons as a major constituent of 
the jet material. In section |i| we discuss synchrotron emis- 
sion, the energy loss processes of the relativistic electrons 
responsible for this emission and a model for the cocoon 
based on KA97 which allows us to calculate its radio emis- 
sion. Section ^| contains a discussion of the results obtained 
with this model for the cocoon and examines changes to the 
evolutionary tracks through the P-D diagram introduced 
by different models for the evolution of the magnetic field in 
the cocoon, the redshift of the source, its orientation with 
respect to the line of sight and the addition of thermal, non- 
radiating particles to the material in the cocoon. 



2 THE RADIO EMISSION OF THE COCOON 

In this section we firstly present a brief review of the syn- 
chrotron emission mechanism and the energy losses affect- 
ing the relativistic electrons responsible for this emission. 
We will assume these electrons to be confined to a small 
volume element. We then go on to develop a model for the 
emissivity of the whole of the cocoon by summing up the 
contribution of each volume element. We neglect the emis- 
sion of the hot spots since we are confining our interest in 
this paper to those radio sources in which the emission of the 
cocoon dominates the total radio luminosity of the source. 



2.1 Synchrotron radiation and loss processes 

The specific volume emissivity due to synchrotron radiation, 
j v , of an ensemble of relativistic electrons in a magnetic field 



of energy density u B averaged over all electron pitch angles 
is (e.g. Shu 1991): 



4 2 2 

-a T cu B /3 e 7 $(v,7)n(7) dj, 



(1) 



where o T is the Thompson cross section, c the speed of 
light, j3 e the speed of the relativistic electrons in units of c 
(/3 e ~ 1), 7 the corresponding Lorentz factor and 71(7) dy the 
number density of electrons with Lorentz factors between 7 
and 7 + 0(7. Taking the average over all possible electron 
pitch angles is justified since we assume that the magnetic 
field is completely tangled on scales much smaller than the 
cocoon itself. The spectral emission of an electron of Lorentz 
factor 7 is given by $(^,7). In this paper we are primarily 
concerned with the total source luminosity, rather then de- 
tails of spectral shape and therefore we use the standard ap- 
proximation that electrons are emitting only at their critical 
frequency v — j 2 u L , where v L is the Larmor frequency; this 
allows us to set $(^,7) = 8(y — ~f 2 v L ). Taking the electron 
ensemble to occupy a volume V the emitted radio power per 
unit frequency and solid angle is 

Pu = ^-o T cu B ^n( 7 ) V. (2) 

If the relativistic electrons are initially accelerated at 
time ti (for example via a first order Fermi process at the 
jet shock) we expect their energy distribution to be a power 
law function of their initial Lorentz factor, 7*, (e.g. Heavens 
& Drury 1988) 



n(~ji,ti) d^i =n ^i i v d'ji 



(3) 



This spectrum will then evolve with time as the electron 
population loses internal energy. The rate of change of the 
Lorentz factor for the electron population is given by 
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The first term on the right hand side represents adiabatic 
expansion losses for the case where the volume V expands 
as V oc t ai (e.g. Longair 1981). The second term is the 
combined loss rate due to synchrotron radiation and inverse 
Compton scattering of the cosmic microwave background 
radiation (CMBR); m e is the electron mass and u c is the 
energy density of the CMBR. Integration of equation (Q) 
gives 



t~ 



7 



7i 
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where for an electron with a Lorentz factor of 7 at time t, the 

Lorentz factor at time ti, when it was initially accelerated, is 

7i, and we may therefore consider 7* to be a function of 7 and 

t. As mentioned above we assume the magnetic field to be 

completely tangled on all relevant scales, allowing us to treat 

u B as a pressure with adiabatic index, F B , such that u B oc 
t -r B a K The 

energy density of the CMBR, u c , is a function 
of redshift (u c oc (z + l) 4 ), but since the lifetime of radio 
sources does not exceed a few times 10 s years (Alexander & 
Leahy 1987) we take u c to be constant during the evolution 
of a given source. Performing the integral in equation (|B|): 



Ct2(t,ti 



(6) 
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where 



d(6U) = Q d(5ti)-p c (ti)d(5V), 



(11) 



Ct2(t,ti) 
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with as = 1 — ai(r B + 1/3) and 04 = 1 — ai/3. If u e {U) is 
the energy density of the relativistic electrons at time U, we 
find for n (ti) from equation (bh 
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where "fi.min and 74, max are the low and high energy cut-off 
of the initial energy distribution. The electrons are uniformly 
distributed over the volume V. Because of the expansion of 
V we have to set t ai n(7,t)d7 = t^nfaijU) d^i. Thus the 
energy distribution of the relativistic electrons at time t can 
be shown to be 



71(7, t) d7 = n a -^5- J c?7, 
where from equation ([]) 

7*7* 



7i 



(9) 



(10) 



t 3 —02(4, *<) 7 
and n is given by equation (Q). 

2.2 A model for the cocoon 

The overall dynamics of the cocoon, and specifically the evo- 
lution of the cocoon pressure, p c , were considered in detail in 
KA97; in this section we review the essentials of this model 
necessary to calculate the luminosity of the source and to 
determine also how an element of radio emitting plasma 
evolves within the cocoon. 

For the calculation of the radio emission of the co- 
coon we divide its contents into three separate 'fluids' with 
individual energy densities. The first fluid consists of the 
electrons (and possibly positrons); these are described by 
the energy spectrum given in equation (0) and contribute 
an energy density u e with adiabatic index T e . For a com- 
pletely tangled field we have a magnetic 'fluid' with en- 
ergy density u B and adiabatic index, F B = 4/3. Finally, 
we allow for the possibility that protons and/or electrons 
with a thermal spectrum are present in the cocoon with 
an energy density u T and adiabatic index F T . The over- 
all dynamics are governed by the pressure in the cocoon, 
Pc = (r c — l)(« e + u B + Ut)', where the adiabatic index of 
the cocoon as a whole, r c , depends on the relative pressures 
of each component. 

To allow for variations of the energy densities within 
the cocoon we split the cocoon into small volume elements 
SV. The fluid present in the volume SV is injected into the 
cocoon at a time t, over a short time interval Sti, during 
which the internal energy of SV, SU, changes according to 



where Q is the power of the jet. During the interval Sti the 
volume elements expand adiabatically which changes their 
pressure from the hot spot pressure Ph(ti) to that of the co- 
coon p c (ti + Sti) ~ Pc(ti). With this assumption integration 
of equation (lll|) gives 



SU = 



p c (ti)8V(ti) 
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From KA97 we take Ph/Pc ~ 4ii 2 which then yields 
SV(U) = (Fc ~ ■ L) . Qo (4i? 2 )^ su. 

Pc \ti ) 

If the expansion of SV after the interval Sti is also adiabatic 
we finally find 



(13) 



8V(t) = 
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where ai = (4 + /3)/(F c (5 — 0)). Here we have used the result 
from KA97 that p c oc t < ' 4 ~ l3)/<5 'i 3) with (3 being the expo- 
nent in the density distribution of the external atmosphere. 

We define k' as the ratio of the energy densities of the 
thermal particles, u T , to that of the electrons, u e , when they 
are injected into the cocoon at U. Note that this differs from 
the usual way of defining k as the ratio of the energy densi- 
ties of non-radiating particles and relativistic electrons. In 
our case the non-relativistic electrons of the power law en- 
ergy distribution (equation ^|) are already included in u e . 
We also introduce r as the ratio of the energy density of the 
magnetic field, u B , to the sum of the energy densities of the 
particles, u e + u T . With these definitions: 



U e (ti 
U B (ti 



p c (ti) 



(r c -l)(fc' + l)(r + l)' 
rpc(ti) 



(15) 



(T c -l)(r + l) 

By identifying V in equation (|^) with SV we can calcu- 
late the radio emission of this volume element, and obtain 
the total emission of the cocoon by summing the contribu- 
tions of all such elements within the cocoon; from equation 
( p^ ) it is clear that this sum reduces to an integration over 
the injection time, ti: 



Pu 



a T cr 



}„n (4R 2 )' 



Q 67T v (r + 1) 
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d2(t,ti) 7 
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(16) 



For the older parts of the cocoon the energy losses of 
the relativistic electrons can be so severe that 74 — > 00 in 
equation (|l^), i.e. electrons with the correct Lorentz factor, 
7, to produce radiation with frequency v now, at time t, 
would have had 7,: — > 00 at injection time, ti. Volume ele- 
ments for which this is the case no longer contribute to the 
radio emission, hence we define a minimum injection time, 
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tmin, at which the cocoon material is still radiating at fre- 
quency v. The integration limits in equation ( Jl6| ) are then 

from train to t. 

The integral in equation ( |l6| ) is not analytically soluble 
for arbitrary p, and we therefore use a numerical approach, 
the results of which we present in Section ^. 



3 EVOLUTIONARY TRACKS 

To perform the calculations discussed in section |^ an ex- 
plicit form for the electron distribution function at the time 
the electrons are injected into the cocoon is required. It is 
shown in KA97 that the material in the jet is 'cold' im- 
plying an adiabatic index of 5/3. Nevertheless, in the rest 
frame of the jet shock the jet material is moving at relativis- 
ts speeds (0.77c — > 0.87c for the parameters used in KA97), 
which implies an exponent for the energy distribution of 
p — 2.14 (Heavens & Drury 1988). This value is derived us- 
ing a test particle approach which may not be appropriate 
for jet shocks in extragalactic sources. However, the spec- 
tral indices of observed radio sources infer 2 < p < 3 (e.g. 
Alexander & Leahy 1987) and we therefore use p — 2.14 for 
the tracks calculated in this section. Since p > 2, all mo- 
ments of the electron energy distribution up to the second 
moment (related to the synchrotron loss rate) are insensi- 
tive to the choice of upper limit for the energy distribution, 
7i,maa;, and specifically converge even as 7i jOTa:c — > oo; we 
therefore take this limiting value. The main form in which 
energy is transported along the jet is as bulk kinetic energy; 
the energy transported per particle in the jet is therefore 
approximately constant. The acceleration process at the jet 
shock transfers this energy to the fraction of the electron 
population constituting the power law energy distribution; 
a substantial proportion of the material in the cocoon must 
therefore be cold, contributing to the density but not the en- 
ergy density. So as to account for this population we assume 
the power law distribution extends to "/i, m in = 1; however 
we shall for the moment neglect any thermal plasma setting 
ft' = 0. 

From KA97 we have expressions for the cocoon pressure 
p c (assuming cylindrical symmetry of the cocoon) and for 
the linear size as a function of time (D oc t 3 ^ 5 ~^). For the 
initial ratio of the energy densities of the magnetic field and 
of the particles we use r — /4 which is taken from min- 

imum energy arguments (e.g. Burbidge 1956). Other typical 
observational parameters are taken as follows. For the den- 
sity profile of the external atmosphere (p x — p (d/a )~^) we 
take the central density p = 7.2 ■ 10 -22 kg/m 3 , for a core 
radius a = 2 kpc and /3 = 1.9; these values are typical for 
a galaxy out to about 100 kpc from its center (Canizares et 
al. 1987). Rawlings & Saunders (1991) find jet powers, Q rs , 
from 10 37 W to 10 39 W for FRII sources. These authors used 
in their calculations a form of the minimum energy argument 
with cut-offs for the power law energy distribution corre- 
sponding to the observational limits of the radio spectrum 
of 10 MHz and 100 GHz respectively implying that there is 
no thermal material present in the cocoon (k — 0). Their 
frequency cut-offs translate in our model to 7i, m i n ~ 500 
and 7i ; max 10 s for a source with D = lOOfcpc. Since in 
our model the power law of the energy distribution extends 
to 7i, m in = 1 there is some thermal material present in the 



cocoon and for the quoted values we find k = 3.6. We must 
therefore adjust their values to be consistent with our choice 
of 7i,min and yi,raax- Rawlings & Saunders (1991) also as- 
sumed that half of the energy that is transported by the jet 
during its life time, t, is lost to the surrounding IGM because 
of the expansion work of the cocoon. However from KA97 
we find that this work is given by 



p c dV c 



5-/3 



9 [r c + (r c - i)i? 2 ] -4-/3 



lot, 



(17) 



where V c is the volume of the cocoon. We can therefore 
derive jet kinetic powers from the values given by Rawlings 
& Saunders (1991) using 



Qo 



9 r c + (r c - i)i? 2 



10-2/3 



P 4 
(ft + 1) 5 +P 



13 Q r 



(18) 



where R — 2, ft = 3.6 and T c = 5/3. 

The pressure in the cocoon also depends on its axial 
ratio, R, and the bulk velocity in the jet, Vj. We adopt R = 
2 which is an average value (Leahy & Williams 1984) and 
Vj = 0.87c implying a Lorentz factor of the bulk motion of 
7j =2. This velocity is a rather arbitrary choice since there 
are no reliable estimates of Vj in the literature. 

The distribution of sources in the P-D diagram is con- 
structed from the data presented by Laing et al. (1983), 
Subrahmanyan et al. (1996) and Cotter et al. (1996); these 
authors quote luminosities at or close to an observing fre- 
quency of 178 MHz and we therefore use this frequency 
in our calculations. Note that the distribution of observed 
sources presented here cannot be taken as the 'true' distribu- 
tion because of the involved selection effects. The observed 
linear size of a radio source projected onto the plane of the 
sky, D, depends of course not only on the advance speed of 
the hot spots but also on sin a, where a is the angle between 
the jet axis and the line of sight. Assuming that the low- 
frequency LRL sample represents a uniform distribution in 
a we adopt the average viewing angle a — 39.5° for the 
model sources. 



3.1 Evolution of the magnetic field 

To proceed we need to know the equation of state of the 
material in the cocoon as a whole and also how the mag- 
netic field is behaving during the expansion of the volume 
elements 8V . Since there is a mixture of three different fluids 
in the cocoon this is not straight forward, and it is also likely 
that the equation of state varies both within the cocoon and 
with time. Hence we investigate three limiting cases: 

• Case 1: Both, the cocoon and the magnetic field energy 
density, have a relativistic equation of state (r c = r B = 
4/3). ' 

• Case 2: The cocoon is 'cold' (F c = 5/3) but the energy 
density of the magnetic field is proportional to the one of 
the relativistic particles and therefore F B = 4/3. 

• Case 3: Cocoon and magnetic field are 'cold' (r c = 
T B = 5/3). This implies some energy dissipation process be- 
tween the magnetic field and the non-relativistic particles 
by which the adiabatic index of the energy density of the 
magnetic field is held constant at 5/3. This effectively keeps 
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Figure 1. Evolutionary tracks for three FRII sources. The upper 
curves are for Q = 1.3 • 10 40 W and 2 = 2, the curves in the 
center for Q = 1.3 • 10 39 W and z = 0.5 and the lower curves 
for Q = 1.3 • 10 38 W and z = 0.2. The other parameters for all 
sets of curves are: R = 2, p = 7.2 ■ 10~ 22 kg/m 3 , a = 2 kpc, 
/3 = 1.9 and k' = 0. Solid curves for case 1, dashed curves for 
case 2 and dotted curves for case 3 (see text). The vertical lines 
on the solid curves are time markers. The first marker on the left 
indicates a life time of 10 6 years of the respective source model. 
Next ten markers for every 10 7 years up to 10 s years. At the last 
marker the source is 2 ■ 10 s years old. Crosses are observed sources 
above 10 25 W Hz -1 sr^ 1 taken from Laing et al. (1983), stars 
from Subrahmanyan et al. (1996) and circles from Cotter et al. 
(1996). 



the energy of the magnetic field and that of the particles in 
equipartition. 

In cases 1 and 3 the ratio of the energy densities of 
the magnetic field and of the particles for every volume ele- 
ment SV in the cocoon, r = 0.785 for the value of p = 2.14 
adopted here, is not changing with time. In case 2 the value 
of r will increase with time since the adiabatic index of the 
magnetic field is lower than the one of the particles. This 
implies that after some time the energy density of the mag- 
netic field will start to dominate the total energy density in 
SV and therefore change the adiabatic index of this part of 
the cocoon. For a volume element injected into the cocoon 
at time U the ratio of energy densities r becomes equal to 1 
at time t ~ 1.9 ■ U. At the same time the total volume of the 
cocoon, V c , has increased by a factor 6.4 (V c oc f 9 /' 5- ' 3 ^ see 
KA97) implying that in most parts of the cocoon r < 1. We 
therefore make the approximation that the adiabatic index 
of the cocoon as a whole is 5/3. 

Figure [l] shows the evolutionary tracks for all three cases 
(corrected for T c = 5/3: Q a = 13 ■ Q ra ) for two limiting jet 
powers Q = 1.3 • 10 38 W and Q = 1.3 • 10 40 W and an 
intermediate case with Q = 1.3 • 10 39 W; the redshift z is 
0.2 for the low power jet, 0.5 for the 'average' jet and 2 for 
the high power jet. 

The range of luminosities observed in FRII sources is 
bounded by the tracks for the two limiting jet powers and 




Length of radio source in kpc 

Figure 2. Comparison of the different loss processes. All model 
parameters as for the intermediate case in figure [l]. The energy loss 
processes of the relativistic electrons considered are different for 
the tracks: Solid curve: adiabatic losses plus synchrotron losses, 
dashed: adiabatic losses plus inverse Compton losses, dot-dashed: 
no losses. 



the drop in luminosity caused by the catastrophic energy 
losses of the relativistic electrons occurs at approximately 
the correct linear size. The importance of radiative losses 
in reducing the source luminosity at large linear sizes (even 
at a low observing frequency) is very significant — the lack 
of large (Mpc) sources at high luminosity is therefore a di- 
rect result of their intrinsic luminosity evolution, and is a 
strong function of frequency. Most of the observed sources 
'crowd' in the region of the diagram between 100 kpc and 
1 Mpc. From the time markers on the tracks in figure [l] 
it is clear that the model discussed here predicts exactly 
this behaviour since model sources spend the longest part of 
their life time in this region of the diagram. The difference 
between cases 2 and 3 is almost negligible. The higher lumi- 
nosity for case 1 results from the fact that in this case less 
of the energy transported by the jet is lost in the expansion 
of the cocoon. Using the appropriate correction factor for 
Qo for this case by setting T c = 4/3 in equation (^), which 
yields Q — 2.3 • Q rs , puts this track very close to the other 
two. Thus we can conclude that the different possibilities for 
the equation of state of the cocoon and for the evolution of 
the magnetic field do not have a significant influence on the 
evolutionary tracks. We will therefore use only case 3 in the 
following. 

Figure ^ shows a comparison of the effects of the var- 
ious loss processes. The most luminous track in this figure 
represents a model where the energy losses of the relativistic 
electrons are completely neglected, i.e. the right hand side 
of equation (^) is set to zero. This curve is equivalent to 
the results of Baldwin (1982). The two other tracks shown 
both incorporate adiabatic losses plus one of the radiative 
loss processes each. Obviously the adiabatic losses lead to a 
pure luminosity off-set which one would expect from equa- 
tion (^) if the two radiative terms on the right are set to zero. 
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1 10 100 1000 10* 
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Figure 3. The influence of the redshift z on the evolutionary 
tracks at v = 178 MHz. Solid curve: z = 0, dashed: z = 0.5, 
dot-dashed: 2 = 1, dotted: z = 2 and dot-dot-dot-dashed: z=3. 
All other model parameters as in the intermediate case in figure 
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Figure 4. The influence of the shape of the density distribution 
in the IGM. Solid curve: (3 = 0, dashed: ,3 = 1, dot— dashed: 
(3 = 1.5 and dotted: (3 = 2. All other model parameters as in the 
intermediate case in figure [l| 



From the figure it is also clear that synchrotron losses are 
dominating the shape of the track for small sources and in- 
verse Compton losses are responsible for the steep decline of 
the luminosity of large sources. The decreasing importance 
of synchrotron losses with increasing linear size is caused 
by the decreasing energy density of the magnetic field in 
the cocoon. Once the energy density of the CMBR becomes 
comparable to the energy density of the magnetic field the 
inverse Compton losses determine the shape of the track and 
limit the size of the radio source. 

As the source luminosity decreases as the result of in- 
verse Compton losses, so the prominence of the hot spots 
will increase. A detailed investigation of this effect requires 
a proper model for the hot spot itself and will be investigated 
in a forthcoming paper. 



3.2 Effects of redshift and environment 

In figure ^| we show evolutionary tracks at an observing fre- 
quency of 178 MHz at various redshifts; the other source 
parameters are given in the figure caption. These tracks are 
calculated assuming that there is no evolution of the sur- 
rounding gas with redshift. Since the energy density of the 
CMBR increases with redshift the inverse Compton losses 
become important earlier and the evolutionary tracks begin 
to steepen for smaller linear sizes of the source. 

The environments of radio sources have a great influ- 
ence on their appearance. The effects on the radio luminosity 
of a source caused by a change in the shape of the density 
profile is shown in Figure ^, where the exponent, (3, in the 
density distribution of the IGM is varied. The track for (3 — 
is purely hypothetical since KA97 have pointed out that a 
jet with the parameters used here could not reach a linear 
size greater than about 80 kpc in an uniform environment 



without being destroyed by turbulence; for this value of (3 
the track is increasing in radio luminosity for small sources, 
a result also found by Baldwin (1982), but the inverse Comp- 
ton losses are still strong enough to reverse this trend and 
cause the track to steepen for large linear sizes. The be- 
haviour for f3 = 1 is close to that of a source with constant 
radio luminosity in the absence of radiative losses. 

The point of intersection is due to the intrinsic char- 
acteristic time scale of the problem r = (a„ po/Qo) 1 ^ 3 (see 
KA97) which is the same for all values of (3. Note also that 
the point at which inverse Compton losses become impor- 
tant shifts to smaller linear sizes for smaller values of (3 (with 
the exception of (3 — 0). This is because the pressure in the 
cocoon is lower at any given linear size in the case of large 
values of (3 compared to smaller values. The population of 
relativistic electrons in these cases has therefore lost less en- 
ergy when the energy density of the CMBR becomes com- 
parable to the energy density of the magnetic field in the 
cocoon. The steepening of the track due to inverse Comp- 
ton losses then occurs at larger linear sizes. In the case of 
/3 = the luminosity of the source would increase with lin- 
ear size for all sizes without any radiative losses. However, 
the synchrotron losses of the electrons in this case are so 
severe that the luminosity of the source is almost entirely 
produced by material which has just been injected into the 
cocoon while older parts are not radiating anymore, lead- 
ing to almost constant radio luminosity for all linear sizes. 
Because of this the exact point at which inverse Compton 
losses become important in this case is difficult to determine 
and the start of the steepening of the track is less obvious. 

The observed distribution of sources in the P-D dia- 
gram seems to be fitted best by tracks with (3 just less than 
2. This is in agreement with (3 = 1.9 found by Cotter (1996), 
who constructed evolutionary tracks through the P-D di- 
agram by comparing the positions of sources with similar 
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Evolutionary tracks of FRII sources through the P-D diagram 7 
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Length of radio source in kpc 

Figure 5. The influence of the central density of the density 
distribution of the IGM. Solid curve: p a = 7.2 ■ 10~ 20 kg/m 3 , 
dashed: p a = 7.2 ■ 10~ 21 kg/m 3 , dot-dashed: p a = 7.2 ■ 1CT 22 
kg/m 3 and dotted: p = 7.2 ■ 10 -23 kg/m 3 . All other parameters 
as in the intermediate case in figure ll 
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Figure 6. The influence of the aspect ratio of the cocoon R. 
Tracks are plotted for Q a = 1.3 ■ 10 39 W, p a = 7.2 • 1CT 22 kg/m 3 , 



a = 2 kpc, /3 = 1.9, k = and z = 0.5. Solid curve: R 
dashed: R = 2, dot— dashed: R = 4 and dotted: R = 6. 



1.3, 



observed jet power, and also with X-ray observations of the 
gas surrounding galaxies at low redshift (Canizares et al. 
1987). 

Figure ^| shows the effect of varying the central density 
of the IGM, p . For a greater central density the pressure 
in the cocoon is higher which leads to an increased radio 
luminosity and also to greater energy losses of the relativistic 
electrons via synchrotron radiation. This effect 'flattens' the 
tracks of small sources in a high density environment. Again 
the steepening of the tracks of large sources occurs at smaller 
linear sizes for sources in denser environments since their hot 
spot advance speeds are smaller. 



3.3 Effects of jet properties 



In section 3.1 we have already investigated the influence of 
the jet power, Q , on the evolutionary tracks. Other jet pa- 
rameters which have an effect on the radio luminosity of 
radio sources are the jet half-opening angle, 9, (which con- 
trols the aspect ratio of the cocoon R) and the nature of the 
jet material. 

If the expansion of the cocoon perpendicular to the jet 
axis is confined by the ram pressure of the IGM we find that 
the aspect ratio of the cocoon is determined by the ratio of 
the pressure at the hot spot to the pressure in the cocoon. 
From KA97 we find that this ratio is proportional to the 
inverse square of the jet half opening angle 9. Figure ^ shows 
the evolutionary tracks for the range of R observed in sources 
at constant Q . The limiting cases R = 1.3 and R = 6 
correspond to 9 — 47.8° and 9 — 10.4° respectively. The 
pressure in the cocoon must be higher in the sources with 
wider cocoons to allow the faster expansion perpendicular to 
the jet axis, leading to higher synchrotron losses and hence 
giving relatively flat tracks for small linear sizes. Since the 




"~ 1 10 100 1000 10 

Length of radio source in kpc 

Figure 7. The influence of protons in the jet. Tracks are plotted 
for p = 7.2- 1CT 22 kg/m 3 , a a = 2 kpc, /3 = 1.9 , z = 2 and R = 2. 
Solid curve: Q a = 1.3 • 10 40 W and k' = 0, dashed: Q = 1.3 10 40 
W and k' = 100 and dot-dashed: Q a = 10 42 W and k' = 100. 



energy density at the hot spots is lower relative to that of the 
cocoons, the hot spot advance speeds are smaller than those 
for the sources with thinner cocoons and the steepening of 
the tracks occurs at smaller linear sizes. 

All models considered up to this point have assumed 
that the jets consist of an electron-positron plasma and that 
the contribution of the thermal particles to the energy den- 
sity in the cocoon is negligible. If there are protons in the 
jet they are also accelerated to relativistic velocities at the 
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jet shock. Bell (1978) showed that for an proton-electron 
plasma accelerated in a shock front the protons can store 
ten times more energy than the electrons. Other accelera- 
tion scenarios predict much higher values (Eilek & Hughes 
1991). From figure | it is clear that an addition of protons 
to the jet material, even if they contribute only moderately 
to the energy density in the cocoon (k' = 100), decreases 
the luminosity of radio sources severely. Should the jets in 
FRII sources consist of protons and electrons they must have 
much higher energy transport rates then we have assumed 
so far in order to explain the most luminous sources. This 
would considerably increase the hot spot advance speeds of 
these sources and lead in turn to very large linear sizes for 
the sources before inverse Compton losses become significant 
(see figure [?]). Since we do not observe such large luminous 
sources this can be taken as evidence against the presence 
of protons in the jet. 



4 CONCLUSIONS 

A model for the cocoon of FRII radio sources based on the 
model of KA97 is presented. Energy loss processes for the 
relativistic electrons producing the radio emission of the co- 
coon are incorporated into a calculation of the expected ra- 
dio luminosity. The resulting evolutionary tracks through 
the P-D diagram are shown to be in good agreement with 
the observed distribution of sources. The lack of luminous 
giants in the diagram is reproduced just by the energy losses 
of the electrons without invoking changes in the density dis- 
tribution of the IGM or the switching off of the central en- 
gine. The exact details of the evolution of the magnetic field 
in the cocoon are shown to have no significant effect on the 
evolutionary tracks. The effects of the environment and jet 
parameters on the tracks are investigated. The requirement 
of very high jet powers and the associated high hot spot 
advance speeds for proton-electron jets to account for the 
most luminous observed sources might rule out large frac- 
tions of protons in the jet. 
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